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ABSTRACT 

As a promising channel to Type la supernovae (SNe la) , we have proposed a symbiotic binary system 
consisting of a white dwarf (WD) and a low mass red-giant (RG), where strong winds from the accreting 
WD play a key role to increase the WD mass to the Chandrasekhar mass limit. However, the occur- 
rence frequency of SNe la through this channel has been still controversial. Here we propose two new 
evolutionary processes which make the symbiotic channel to SNe la much wider. (1) We first show that 
the WD -f RG close binary can form from a wide binary even with such a large initial separation as 
a-i < 4OOOOi?0. Such a binary consists of an AGB star and a low mass main-sequence (MS) star, where 
the AGB star is undergoing superwind before becoming a WD. If the superwind at the end of AGB 
evolution is as fast as or slower than the orbital velocity, the wind outflowing from the system takes 
away the orbital angular momentum effectively. As a result the wide binary shrinks greatly to become a 
close binary. Then the AGB star undergoes a common envelope (CE) evolution. After the CE evolution, 
the binary becomes a pair of a carbon-oxygen WD and the MS star. When the MS star evolves to a RG, 
a WD + RG system is formed. Therefore, the WD + RG binary can form from much wider binaries 
than our earlier estimate which is constrained by Ui <^ ISOOi?©. (2) When the RG fills its inner critical 
Roche lobe, the WD undergoes rapid mass accretion and blows a strong optically thick wind. Our earlier 
analysis has shown that the mass transfer is stabilized by this wind only when the mass ratio of RG/WD 
is smaller than 1.15. Our new finding is that the WD wind can strip mass from the RG envelope, which 
could be efficient enough to stabilize the mass transfer even if the RG/WD mass ratio exceeds 1.15. If 
this mass-stripping effect is strong enough, though its efficiency is subject to uncertainties, the symbiotic 
channel can produce SNe la for a much (ten times or more) wider range of the binary parameters than 
our earlier estimation. With the above two new effects (1) and (2), the symbiotic channel can account for 
the inferred rate of SNe la in our Galaxy. The immediate progenitor binaries in this symbiotic channel 
to SNe la may be observed as symbiotic stars, luminous supersoft X-ray sources, or recurrent novae like 
T CrB or RS Oph, depending on the wind status. 

Subject headings: binaries: symbiotic — stars: individual (T CrB, RS Oph) — stars: mass-loss — 
stars: novae — supernovae: general — X-rays: stars 



1. INTRODUCTION 



It is widely accepted that Type la supernovae (SNe la) 
arc thermonuclear cxploaiona of accreting white dwarfa 



(e.g., Nomoto, Iwamoto, fc Kiehimoto 1007). — However, 
whether the explosion of the white dwarf takes at the 
i Chandraackhar maaa limit or at the aub Chandraackhar 



mass has been controversial (e.g., Nomoto et al. 1004; Ar 



nett 1996; Branch 1998). Also the issue of double degen 
erate (DD) vs. single degenerate (SD ) for the progenitor 
scenario has been still debated (e.g.. Branch et al. 1995 
for a review). 

The DD scenario assumes that merging of double C+0 
white dwarfs with a combined mass surp assing the Chan- 
drasekhar mass limit induces SN la (e.g., Iben & Tutukov 
1984; Webbink 1984). However, this scenario has not been 
well supported. Observationally, the search for DDs has 



discovered only several systems whose combined mass is 
less than the Chandrasekhar mass or wh ose separation 



is too wide to merge in a H ubble time ( Branch et al 
1995; Renzini 1996; Livio 199^ for reviews). Theoretically, 



the DD has been suggested to lead to accretion-induced 
collapse rathe r than SN la (Nomoto fc Iben 1985; ^aio & 



Nomoto 1985|, 1998; [Segretain et al. 19971) 



For the SD scenari o, observational counterp arts may be 
symbiotic stars (e.g., Munari & Renzini 1992). Kenyon et 
al. (1993) and Renzini (1996) have suggested that sym- 
biotics are more likely to lead to the sub-Chandrasekhar 
mass explosion because the available mass in transfer 
may not be enough for white dwarfs to reach the Chan- 
drasekhar mass. However, photometric and spectroscopic 
features of majority of SNe la are in much better agree- 
ment with the Chandrasekhar mass model than the sub- 
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Chandrasekhar ma ss model ( Hoflich fc Khokhlov 1996 ; 
Nugent et al. 1997| ). 

Hachisu, Kato, & Nomoto (1996; hereafter HKN96) 
have shed new hght on the above SD/sub-Chandrasekhar 
mass scenario and proposed a new progenitor system based 
on the optically thick wind theory of mass-accreting white 
dwarfs. In the scenario of Iben & Tutukov (1984) and Web- 
bink (1984), they excluded close binary systems consisting 
of a mass-accreting white dwarf (WD) and a lobe-filling 
red(sub)-giant (RG) (first discussed by Whelan fc Iben| 
1973|] mainly because such a system suffers from unstable 
mass transfer when the mass ratio of the RG/WD exceeds 
0.79, i.e., q = Mrq/Mwd > 0.79. However, HKN96 have 
shown that optically thick winds from the mass-accreting 
white dwarf stabilize the mass transfer up to g < 1.15 even 
if the RG has a deep convective envelope. Such an object 
may be observed as a symbiotic star. Thus they proposed 
a new channel to SNe la, through which a white dwarf ac- 
creting mass from a lobe-filling red-giant can grow up to 
the Chandrasekhar mass limit and explodes as an SN la 
(SD/Chandrasekhar mass scenario of symbiotics). 

Li & van den Heuvel (1997) reanalyzed the HKN96 
model and identified two isolated regions for SN la pro- 
genitors in the initial orbital period vs. the initial donor 
mass plane, i.e., in the logPg ^ -^^d,o plane: 1) One is 
a relatively compact close binary consisting of a Afd.o '^ 
2 — 3Mq slightly evolved main-sequence companion and a 
A^WD.o '^ 1.0 — 1.2AIq white dwarf with the initial orbital 
periods of Pq ~ 0.5 — 5 d (hereafter, WD-fMS systems), 
and 2) the other is a relatively wide binary consisting of 
a low mass (Md,o '^ l-^o) red-giant companion and a 
-^WD.o ~ 1.2Mq white dwarf with the initial orbital pe- 
riods of Po ~ 100 - 800 d (hereafter, WD-hRG systems). 
They also concluded that the new model accounts for the 
inferred rate of SNe la in our Galaxy. 

In Li & van den Heuvel's (1997) analysis, their SN la 
progenitor region for the WD+RG system is very small 
compared with the region for the WD-I-MS system. The 
contribution of the WD-I-RG systems to the total rate of 
SN la explosions was expected to be very small or negli- 
gible. It is because the WD-f-RG channel is restricted by 
the condition for the stable mass transfer, i.e., q < 1.15. 
In the present paper, we propose a new evolutionary pro- 
cess which makes the WD-I-RG channel to SN la much 
wider than that of IIKN96's original modeling. We in- 
clude a mass-stripping effect of the red-giant envelope by 
the wind. The effect removes the limitation oi q < 1.15 
and, as a result, the new parameter region producing an 



three reasons why our estimates are much larger than Yun- 
gelson & Livio's (1998) estimates: 1) We introduce a mass- 
stripping process into the WD-I-RG systems as an exten- 
sion of IIKN96's model. As a result, the new parameter re- 
gion producing an SN la becomes ten times or more wider 
than the previous region calculated by HKN96 and Li & 
van den Heuvel (1997) when the efficiency of the mass- 
stripping effect is strong enough. 2) Yungelson & Livio 
(1998) assumed that the initial separation is smaller than 
Qi ^ 1500 Rq in their estimations of the WD+RG model. 
If one includes the effect of angular momentum loss by slow 
winds at the end of stellar evolution, however, very wide 
binaries with the separation of 1500 Rq <, ai < 40000 Rq 
shrink into 30 Rq <, o,f ^ 800 Rq, which provide ap- 
propriate initial conditions for our WD+RG models. 3) 
We believe that Yungelson & Livio (1998) did not include 
an important evolutionary path in the rate estimation for 
Li & van den Heuvel's (1997) WD+MS model. They as- 
sumed that relatively massive white dwarfs (~ IMq) are 
born from an AGB star and neglected the possibility that 
it comes from a helium star whose hydrogen-rich enve- 
lope is stripped away in a common envelope evolution at 
the red-giant phase with a helium core. Very recently, in- 
cluding the evolutionary path mentioned above, Hachisu, 
Kato, Nomoto, & Umeda (1999, hereafter HKNU99) have 
shown that the realization frequency for the WD+MS sys- 
tems is as large as i^ms ~ 0.001 yr~^, which accounts for 
one third of the inferred rate of SNe la in our Galaxy. 
It has been argued that s ome of the recurrent n ovae are 

because 



progenitors of SNe la (e.g., Starrfield et al. 
these white dwarfs are suggested to be very massive and 
close to the Chandrasekhar mass limit. Morphologically, 
recurrent novae are divided into three groups according to 
the companion star; dwarf companions, slightly evolved 
main-sequen ce (or sub-giant) companion s, and red-giant 
companions (Schaefer & Ringwald 1995). The latter two 
groups are relevant to SNe la progenitors and the good 

examples are as follows. 

1) T CrB (Porb = 22 7.67 d; e.g., [Lines et al. 1988| ) and 
RS Oph (Porb = 460 d; [Dobrzycka fc Kenyon 19l1| Tbelong 
to the last group of red-giant companions and their white 
dwarf mass es are very close to the Chandrasekhar mass 
limft (e.g., iKato 1990. 1995, 1999: IShahbaz et al. 1997; 



Belczyhski fc Mikolajewska 1998| ; [Hachisu fc Kato 1999a|) 



SN la becomes ten times or more wider than the prcvi- 



These two systems correspond to the WD+RG systems 
These extremely massive white dwarfs are naturally ex- 
plained in our SN la progenitor scenario. 

2) On the other ha nd, U Sco (Po,b = 1-23 d; [Schae- 



ous region calculated by HKiN96 and Li & van den Heuvel 
(1997). 

Recently, Yungelson fc Livio (1998) claimed, based on 
their population synthesis results, that HKN96's and Li 
fc van den Heuvel's (1997) model can account for only, at 
most, 10% of the inferred rate of SNe la. Introducing a new 
evolutionary process into HKN96's modeling, we have re- 
analyzed the SN la rate for our extended HKN96's model 
and also for Li & van den Heuvel's model. Our present 
analysis reveals that the realization frequencies of SNe la 
coming from our WD+RG/WD+MS models are ~ 0.002 
yr^^/^ 0.001 yr~^, respectively, and the total SN la rate 
becomes ^ 0.003 yr~^, which is large enough to account 
for the inferred rate of SN la rate in our Galaxy. There are 



fer fc Ringwald 1995[ ) and V394 CrA {Porh = 0.755~a', 
Schaefer 1990| ) belong to the middle group of the slightly 
evolved main-sequence companions. For this group, it 
has been suggested that the companion has an extremely 
helium-rich envelopes and the primary is a very massive 
white dwarfs close to the Chandrasekhar mass limit. Our 
WD+MS model yields the secondary star having a helium- 
rich envelope as well as the primary of very massive white 
dwarfs as suggested first by Hachisu & Kato (1999b). 

In §2, we describe a new idea of mass-stripping effect 
by strong winds and the mass accumulation efficiency. We 
then search for the initial parameter regions that can pro- 
duce SNe la in §3. In §4, we discuss relevance to recurrent 
novae, our criticism to Yungelson & Livio's claims, esti- 
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mation of SN la rates of the WD+RG/WD+MS systems, 
and the possibihty of detecting hydrogen- Hnes at SN la ex- 
plosions during the strong wind phase. Conclusions follow 
in §5. 

2. PROGENITOR SYSTEMS 

First of all, we illustrate a full evolutionary path of our 
WD-fRG system from the zero age main-sequence stag e 



{stage A) to the SN la explosion (stage F) in Figure Al 

A) Zero age main-sequence. 

B) The primary has evolved first to become an AGB 
star and blows a slow wind (or a super wind) at the 
end of stellar evolution. 

C) The slow wind carries the orbital angular momentum 
and, as a reaction, the separation shrinks consider- 
ably (by about a factor of ten or more), which is a 
similar process to the common envelope evolution. 

D) A carbon-oxygen white dwarf (the initial primary) 
and a zero age main-sequence star (the initial sec- 
ondary) remain. 



E) p'hc initial secondary has evolved to a red-giant 
forming a helium core and tills up its inner critical 
Roche lobe. Mass transfer begins. The WD compo- 
nent blows a strong wind and the winds can stabilize 
the mass transfer even if the RG component has a 
deep convective envelope. 

F) The WD component has grown in mass to the Chan- 
drasekhar mass limit and explodes as a Type la su- 
pernova. 

For an immediate progenitor system of Type la su- 
pernovae (SNe la), we consider a close binary initially 
consisting of a carbon-oxygen white dwarf (C-l-0 WD) 
with MwD.o = 0.6 — 1.2Mq and a low-mass red-giant 
star with A/rq.o = 0.7 — S.OM© having a helium core 
of AfHo.o = 0.2 - 0.46AfQ {stage E). The initial state of 
these immediate progenitors is specified by three param- 
eters, i.e., MwD.O: -^rg.Oj and the initial orbital period 
Pq {M]ie,o is determined if Pq is given). We follow binary 
evolutions of these sy stems by using empirical formulae 
(Webbink ct al. 1982) and obtain the parameter range(s) 
which can produce an SN la. 

2.1. Conventional evolution scheme 

When the companion evolves to a red-giant (RG) and 
fills its inner critical Roche lobe, mass transfer begins from 
the RG to the WD. If both the total mass and the total 
angular momentum are conserved and the mass transfer is 
steady, its rate is given by 



M2~\r2 



/H{q), 



(1) 



EV 



where {R2/ R2)wv represents specifically the evolutionary 
change in the secondary radius and 

^(.)-^(l + .)-2(l-,), (2) 

ding 



where q is the mass ratio defined by 

gEEMs/Mi, 



(3) 



{All is the mass of the primary, i.e., the WD component, 
and M2 the mass of the secondary, i.e., the RG compo- 
nent). Here we use the empirical formula proposed by 
Eggleton (1983), 



R% 



fiq) 



0.49fl2/3 



0.6g2/3 + ln(l + qi/3) 



(4) 



for an effective radius {R2) of the secondary's inner critical 
Roche lobe. For the separation a, we simply assume a cir- 
cular orbit. To estimate (P2/P2)ev we use the empirical 
formulae proposed by Webbink et al. (1983). 

For a sufficiently large mass of the secondary M2 (i.e., 
q > 0.79), however, equation (1) gives a positive value of 
M2 ■ This means that the mass transfer proceeds not on an 
evolutionary time scale but rather on a thermal or dynam- 
ical time scale. The gas falls very rapidly onto the WD and 
forms an extended enve lope around the WD (e.g 
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et al. 1979| ; [ben 19881 ). This envelope expands to fill the 
inner and then outer critical Roche lobe. It eventually re- 
sults in the formation of a common envelope, in which the 
two cores are spiraling in each other. It forms a very com- 
pact binary system consisting of a C-l-0 WD and a helium 
WD, or a merger of C 4-0 and He cores. These systems 
have been extensiv ely examined by many authors (e.g., 
Iben & Livio 1993 for a review and references therein). 



2.2. White Dwarf Winds 



Ho wever, the recent version of opacity ( Iglesias fc Rogers 
1996) has changed the story. Optically thick winds are 



driven when the WD envelope expands and the p hoto- 
spheric temperat ure decreases below logTph ~ 5.5 (Kate 
fc Hachisu 1994 ). We have calculated such wind solu- 
tions for various white dwarf masses of Mwd — 0.6, 0.7, 
0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.35 and 1.377Mq and show 
six of ten cases, i.e, Mwr> = 0.6, 0.8, 1.0, 1.2, 1.3 and 



1.377Mq in Figures A2-A6. Here, we choose 1.377M0 as 



a limiting mass just be low the mass at the SN la explo- 
sion in W7 (1.378Mq, [Nomoto et al. 1984[ ) as was done 
in Kato (1995, 1999) 
opacity ( Iglesias 
near logT 



We have used the updated OPAL 

Rogers 1996| ) because its strong peak 

5.2 is about 20 — 30% larger than that of 
the original OPAL opacity ( Rogers fc Iglesias 1992| ) which 
was used in HKN96. The numerical method and various 
assumptions are the same as in Kato & Hachisu (1994) so 
that we omit the details of the numerical calculations to 
avoid the duplicati on. The only difference between the re- 
sults in HKN96 (or Kato & Hachisu 1994) and the present 
ones is the opacity as mentioned above. 

Each wind solution is a unique function of the envelope 
mass AM if the white dwarf mass is given. The envelope 
mass is decreasing due to wind mass loss Af-„,ind (< 0) and 
hydrogen shell burning Mnuc (< 0), i.e.. 



—AM = -M2 
dt 



M 



wind 



M„ 



(5) 
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When the mass transfer rate from the companion il/2 (< 0) 
does not vary so much in a thermal time scale of the WD 
envelope, the WD envelope reaches an equilibrium of 



M2 = Mwind + Mn 



(6) 



Thus we regard the ordinates in Figures A2-A6 as the 
mass tra nsfer rate from the companion. 



Figure A2 shows the envelope mass of the wind and 



static solutions against the total mass decreasing rate of 
the envelope \dM /dt\, i.e., the mass transfer rate from the 
companion \M2\- There exists only a static (no wind) 
solution below the break of each solid line while there 
exists only a wind solution above the break for a given 
envelope mass of AM. The optically thick winds blow 
when the mass transfer rate from the companion star 

IM2I = |M„uc + Mwind I exceeds 



M„ w 0.75 X 10"^ ( -77-^ - 0.40 j Mq yr" 



(7) 



which is reduced fro m ou r wind/static solutions as shown 
in Figure A2 . Figure A3 shows the photospheric tempera- 
ture against the mass transfer rate. It should be noted 
that the optically thick winds begin near logT ^ 5.5, 
which corresponds to the shoulder of the strong peak of 
OPAL opacity in the high temperature side. Figure A4 
depicts the photospheric radius against the mass transfer 
rate. The optically thick winds occur when the photo- 
sphere expands to i?ph ~ O.li?©, which is much smaller 
than the inner critical Roche lobe. We plot the photo- 
spheric velocity v-p\i in Figure ^^ and the ratio Wph/^csc 
between the photospheric vel ocit y and the escape velocity 
at the photosphere in Figure A6 against the mass transfer 
rate. Here, we call the wind as "strong" when the photo- 
spheric velocity exceeds the escape velocity there. When 
the wind is strong enough, the photospheric velocity is as 
high as ~ 1000 km s^^ being much faster than the orbital 
velocity, i.e., Wph > afiorb- 

The optically thick wind is a continuum-radiation driven 
wind in which the acceleration occurs deep inside the pho- 
tosphere (e.g., Kato & Hachisu 1994). Further accelera- 
tion of the wind near the photosphere is negligibly small 
because the photon momentum near the photosphere is 
much smaller than the momentum of the wind. Therefore, 
our results are almost independent of the treatment of the 
radiative transfer near the photosphere. In other words, 
the critical points of the wind solutions, where the wind 
is accelerated, exist deep inside the photosphere and also 
deep inside the inner critical Roche lobe for our WD+RG 
systems. This means that the binary motion does not af- 
fect the acceleration of the optically thick winds mainly 
because the strong winds are already accelerated deep in- 
side the inner critical Roche lobe up to the velocity much 
faster than the orbital motion. 

Once the strong wind occurs, the mass transfer rate is 
modified as 



t).--Ht 




Hi{q) = - 



dlnfjq) 
dlnq 



1 



i + q 



- 2 + 2£ 




H2{q) = 



d\nf{q) 
dlnq 



1 + g 



-2 + 2C(l + 9), 



(10) 



where £^ is the specific angular momentum of the wind in 
units of a^riorb, 



(11) 



J the total angular momentum, M the total mass of the 
system, and florh the orbital angular velocity. The wind 
velocity is about several hundred to one thousand km s~^ 
for relatively massive WDs when the mass transfer rate 




A5). The wind velocity is 



IS ^ 1 X 10-^Mq yr-i (Fig. 

about ten times faster than the orbital velocity because of 
al^orb -- 30 - 100 km s"! for AfwD ~ IMq, M2 ~ IMq, 
and a ~ 30 — 400 Rq. In such cases, winds cannot get an- 
gular momentum from the orbital motion by torque during 
its journey, so that wind has the same specific angular mo- 
mentum as the WD, which is estimated as 



/ — 



1 



(12) 



In this case, function H2{q) changes its sign at q — 1.15. 
Wind mass loss stabilizes the mass transfer in the region 
from q = 0.79 to q = 1.15 (HKN96). 

2.3. Growth of white dwarfs 

During the strong wind phase, as hydrogen steadily 
burns on the surface of the WD, the WD accretes the 
processed matter approximately at a rate of Ma, which 
is given in equation (Mj. When the mass transfer rate de- 
creases below this critical value, optically thick winds stop. 
If the mass transfer rate further decreases below 



M«< 



1 



:Afcr, 



(13) 



which is also reduced from our solutions for the minimum 
envelope mass cor responding to the lower end of each solid 
lines in Figure A2 , hydrogen shell burning becomes unsta- 
ble to trigger very weak shell flashes. Once a shell flash 
occurs, a part of the transferred hydrogen-rich matter may 
be blown off and we need to estimate the net mass accu- 
mulation in order to examine whether or not the WD will 



grow up to th e Chandrasekhar mass limit (e.g., Kovetz & 



Prialnik 1994). Here, we roughly assume that all the pro- 



cessed matter is retained until the accretion rate becomes 

below 

1 



A'f, 



low 



= -A/c 
8 



1 X 10" 'Mq yr" 



(14) 



This assumption may slightly underestimate the mass loss 
from the WD, but there still exists a large uncertainty in 
relation to the mass accumulation ratio, 7/h, i-e., the ratio 
of the retained mass to the transferred mass after many 
cycles of hydrogen shell flashes. To summarize, we assume 
the mass accumulation ratio of hydrogen shell burning as 



0, for|M2|<Miow 
r]H = <; 1, for Mow < IM2I < M„ 

1-Mw/M2, forMer<|M2| 



(15) 



^ Mhig 
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where M„ is the wind mass loss rate calculated by the opti- 
cally thick wind theory and Afhigh '^ 1 x 10^'* — 1 x 10~^Mq 
yr~^ is the upper limit of the ma ss t ransfer rate for our 
wind solutions as shown in Figure AS . 

The steady hydrogen shell burning converts hydrogen 
into helium atop the C+O core and increases the mass 
of the helium layer gradually. When its mass reaches a 
certain value, helium ignites. For the accretion rate given 
by equation (|^), helium shell burning is unstable to grow 
to a weak flash. Once a helium shell flash occurs on rel- 
atively massive white dwarfs (Mwd ^ 1-2M0), a part of 



the envelope mass is blown off in the wind ( Kato et al\ 
1989P] Recently, Kato & Hachisu (1999) have recalculated 
the mass accumulatio n efficiency for helium s hell flashes 
with the new opacity ( Iglesias fc Rogers f996| ). Here, we 
adopt their new results in a simple analytic form, 



??Hg 



-O.fTsflogAfHc + S.Ss") -|-f.05, 

for -7.3 <logMHc < -5.9 
1, for -5.9<logA/Hc^-5 



(16) 



where the mass accretion rate, A/hc, is in units of Mq yr ^ 
and given by 

Mhc-??h-|M2|, (17) 

when hydrogen-rich matter is transferred from the com- 
panion star. Eventually, we have the growth rate of C+O 
WDs given by 



A/c+o = ?/Ho • A/ho = mc ■ m ■ \M2\, (18) 

in our WD-I-RG systems. We use formula (Iq) for various 
white dwarf masses and accretion rates, although their re- 
sults are give n only for a 1.3Mq white dwarf (Kato & 
Hachi fju 1999| ). 

The wind velocity in helium shell flashes reaches as high 
as ^ 1000 km s^^, which is much faster than the orbital ve- 
locities of our WD+MS binary systems, i.e., afloih ^ 300 
km s~^ or of our WD-I-RG binary systems, i.e., aQoih ~ 30 



km 



ft should be noted here that either a Roche lobe 



overflow or a common envelope does not play a role as a 
mass ejection mechanism because the envelope matter goes 
away quickly from the system without int eracting with the 
orbital motion (see Kato & Hachisu 1999 for more details). 



2.4. Mass-stripping effect 

We propose here a new effect, namely, stripping of a 
red-giant envelope by the wind, which is not included in 
HKN96's modeling. Very fast strong winds collide with 
the surface of the companion as illustrate d in Figure Al 
{stage E) and, for more detail, in Figure A7. The red- 
giant surface is shock-heated and ablated in the wind. We 
estimate the shock-heating by assuming that the velocity 
component normal to the red-giant surface is dissipated 
by shock and its kinetic energy is converted into the ther- 
mal ener gy o f the surface layer of the red-giant envelope 
(see Fig. A7). The very surface layer of the envelope ex- 
pands to be easily ablated in the wind. To obtain the mass 
stripping rate, we equate the rate times the gravitational 
potential at the red-giant surface to the net dissipation 
energy per unit time by the shock, i.e., 

5"^ (dA^2) ^ r]cs ■ -v"^ sin^ (3 ■ pv sin f3 

rl2 ^ ^ strip Z 



X 2TTR^smed9, 
between 9 and 6 + d9. where 



a + /3 + 9 = -, 



i?2 _ a 

sin a sin(7r/2 -f /3) ' 



(19) 

(20) 
(21) 



cos "^ I — cos/3 
a 



(3 



-i(/(<z)-cos/3)-A 



(22) 



under the condition of /3 > (See Fig. A7). Here v is the 
wind velocity, R2/a is replaced with f[q) given in equa- 
tion(H), and rj^s is a numerical factor representing the ef- 
ficiency of the energy which is transferred to ablation. We 
assume rjes ~ 1 in the present calculation but examine the 
case of rjcs = 0.3 for comparison. By integrating over 9 
and using the wind mass loss rate M^ (negative) given by 



Afw = -Airr^pv, 



we have 



GM2 



[M2) 



strip 



= 7:V^M„ ■ 77eff ■ g{q). 



(23) 



(24) 



Here g{q) is the geometrical factor of the red-giant sur- 
face hit by the wind including dissipation effect (Fig. A7) 
and it is only a function of the mass ratio q. Then, the 
stripping rate of the red-giant envelope is estimated as 



M, 



[M2) 



9il) 



strip 

13=0 
/3=7r/2 



Vc« ■ 



2GM' 



■M^-giq), (25) 



1 7?^ 
-sin^ P ^ sin 9d9, 

2 r^ 



(26) 



where some numerical values oi g{q) are given in Table Al 



The total angular momentum loss rate by stripping is also 
estimated as 



(■') 



strip 



?7cff 



V'R2 
2GM2 



Afwa^Oorb • h{q), (27) 



Hq) 



/3=0 

/3=7r/2 
1 



1 /?^ 

-sin^p^sin9d9 

2 r^ 



R2 



cos 9 




^ I = 4a'r!„,K = -^ • a^n^ 



strip 



5(9) 



(28) 



(29) 



where the ablated gas is assumed to have the angular mo- 
mentum at the red-giant surface. We note that the nu- 
merical factor 4 = h{q)/g{q) given in Table Al is rather 
small compared with £^ in equation (f2h. 
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Including the effect of mass-stripping, we estimate the 
mass transfer rate from the secondary to the primary 
(Mt < 0) as 



HM = H,iq) + -H2{q) + 2 (4 - ^w) (1 + g)-, (31) 



Ms as summarized in Table A2. Thus two thirds 



H4{q) = (1 + Ci)H2{q) + 2 (4 - ^w) (1 + q)ci, (32) 

?2 



V'R2 



c\ = ?7off • 



■5(9)> 



(33) 



2GM2 

where c\ is a numerical factor indicating stripping effect 
as 

il, = ciM^, (34) 

and estimated as 



ci ~ ?7cff 



.9(9) 



0.0257 V 1000 km s" 
/ R2 \ f JVh^ "^ 



VSOEq/ V2M(7 



(35) 



The stripping effect is important, i.e., ci ~ 1, when the 
orbital period is longer than P ~ 30 d for rjcs — 1, 



WD 



I.OMq, M; 



2.OM0, and V = 1000 km s" 



Function Hi{q) remains negative even for q > 1.15 because 
the second term in the right hand side of equation (32) is 
always negative, i.e., is < iw Therefore, the mass strip- 
ping effect stabilizes the mass transfer even for q > 1.15, 
so that the limitation of g < 1.15 proposed by HKN96 for 
stable mass transfer is removed. 

3. RESULTS 

Our progenitor system can be specified by three ini- 
tial parameters: the WD mass Mwd.Oj the red-giant mass 
Mrg,o, and the orbital period Pq. We study three cases 
P1-P3 (P stands for Progenitor) of such close binary evolu- 
tions. We start the calculation when the secondary fills its 
inner critical Roche lobe. The initial p aram eters (AfvvD.o? 
Mrg.0 7 -Pq) are summarized in Table A2. For example, 
MwD,o = I.OM0, Mrg,o = 2.OM0, and Pq = 300 d for 
case PI. Th e evolutionary histories are plotted in Fig- 



ures A8- AIO , In these three cases, the WDs grow up to 
Aiia — I.38M0 to trigger an SN la explosion as follows. 

PI) The mass transfer begins at a rate of Alt — —8.7 x 
1O~^M0 yr~^, and then the WD wind starts to blow 
at a mass loss rate of My, = —4.8 x 1O~^M0 yr~^; 
and the wind then induces the mass stripping at a 
rate of Ms = -2.0 x 1O-*'M0 yr"!. Thus a half 
of the transferred matter is blown off in the wind. 
The mass transfer rate gradually decreases because 
of decreasing M2/MWD but still higher than Mcr just 
when the WD mass reaches Afja and explodes as an 
SN la at i = 7.2 x 10^ yr during the WIND phase 



(Fig. A8) 



P2) The mass transfer, the WD wind, and the mass strip- 
ping start as in case PI at a rate of Mt, Mw, and 



of the transferred matter accumulates onto the WD. 
The mass transfer rate is gradually decreasing and 
becomes lower than Mcr at t = 6.7 x 10^ yr. The 
wind stops but nuclear burning is still stable until 
the WD mass reaches Mi a to explode as an SN la at 
i = 7.4 x 10^ yr (Fig. ^. After the wind stops, the 
progenitor may be observed as a luminous super-soft 
X-ray source (SSS). 

P3) The mass transfer, the WD wind, and the mass 



stripping start at the rate summarized in Table A2 



In this case, a large part of the transferred mat- 
ter is accumulating onto the WD. The mass trans- 
fer rate decreases and becomes lower than Mcr at 
t = 5.3 X 10^ yr. The wind stops but the nuclear 
burning is still stable until t — 6.5 x 10^ yr. The pro- 
genitor may be observed as an SSS during this steady 
hydrogen shell burning phase. Then nuclear burning 
becomes unstable to trigger very weak shell flashes 
but most of the processed matter accumulates onto 
the WD. The progenitor may be observed as a re- 
current nova (RN) . The WD mass eventually reaches 
Mia to explode as an SN la at i = 1.04 x 10^ yr (Fig. 



AlO) 



We have three cases of the immediate progenitors of our 
WD+RG systems corresponding to cases P1-P3, which are 
summarized in Table A3, i.e., PI) wind phase (denoted by 
WIND), P2) steady nuclear burning phase (denoted by 
SSS), and P3) unsteady weak shell flash phase (denoted 
by RN). In cases P2 and P3, the SSS phase is rather short 
compared with the wind phase. Therefore, the frequencies 
of the SSS phase may be small in our symbiotic channel 
to SNe la. 

The final outcome of the binary evo lutio n is summa- 
rized in the logPo — -^^d,o plane (Fig. All). The right 
region of long Pq in the figure represents our new re- 
sults of the WD-I-RG system. In the left region of short 
Pq, we also show the results of the WD-I-MS system for 
comparison (see IIKNU99 for details). Each grid in the 
log Pq — Aid.o plane corresponds to the evolutionary model 
of our wide WD+RG systems (labeled by "WD-I-RG sys- 
tem") together with the compact WD-fMS systems (la- 
beled by "WD-hMS system"; HKNU99). Here, Md,o is 
the mass of the donor, i.e., Mrg,o (the initial mass of the 
red-giant component), or Mms.o (the initial mass of the 
slightly evolved main-sequence component). The initial 
mass of the white dwarf is assumed to be AfwD.o = I.OM0. 

The outcome of the evolution at the end of our calcula- 
tions is classified as follows. 

i) Formation of a common envelope where the mass 
transfer is unstable {Hi(q) > 0) at the beginning of 
mass transfer (denoted by x). 

ii) SN la explosions (denoted by ® corresponding to case 
PI, O corresponding to case P2, and correspond- 
ing to case P3), where the WD mass reaches 1.38 

Mq. 

iii) Novae or strong hydrogen shell fiash(denoted by 
open triangle) , where the mass transfer rate becomes 
below Afiow. 



Hachisu, Kato, and Nomoto 



iv) Helium core flash of the red giant component (de- 
noted by filled triangle) where a central helium core 
flash ignites, i.e., the helium core mass of the red- 
giant reaches O.46M0 . 

The region enclosed by the thin solid line produces SNe la. 
In HKN96's model, this region was limited by Mrg.o < 
1.15Mq for MwD.o = 1-Mq. This new area is about ten 
times or more wider than that of HKN96's modeling. 

We also show, in Figure A12| , other three cases of the ini- 
tial WD mass, AfwD,o = O.SMq, 0.9Mq, and I.IMq (thin 
solid) together with Mwd.o = 1-OMq (thick solid). The 
regions of Mwd,o = 0.6Mq and 0.7Mq vanish for both the 
WD-I-MS and WD-I-RG systems. It is clear that the new 
region of the WD-I-RG system is not limited by the condi- 
tion oi q < 1.15, thus being ten times or more wider than 
the region of HKN96's model for the other initial white 
dwarf masses. 

For MwD,o > I.2M0, the central density of the WD 
' 10 -■^ 



reachtfs '^ 10 g cm bcfuie hcaLiug wave fium Llie hy- 
droge n burning layer reaches the center. As a result, the 



-^^WD = 1.2 ± 0.2Mq by Belczyhski and Mikolajewska 
(1998) and the other is Mwd = 1.3 - 2.5Mq by Shah- 
baz et al. (1997). Belczyiiski and Mikolajewska de- 
rived a permitted range of binary parameters from am- 
plitude of the ellipsoidal variability and constraints from 
the orbital solution of M-giants. The white dwarf mass 
is permitted up to 1.44Mq under the condition of a cer- 
tain mass ratio and inclination of the orbit (in their Fig. 
4). In Shahbaz et al. (1997), a massive white dwarf of 
-A^WD — 1-3 — 2.5Mq in T CrB is also suggested from 
the infrared light curve fitting. Combining these two per- 
mitted ranges of the white dwarf mass in T CrB, we 
may conclude that a mass of the white dwarf is between 
A^WD = 1.3 — I.4M0, which is very consistent with the 
hght curve analysis (AfwD ^ 1.37- 1.38A/q of T CrB) by 
Kato (1999). 

However, it is very unlikely that such very massive white 
dwarfs were born at the end o f single star evolut ion in a 
binary (e.g., Weidemann 1986; see also eq.(2) of Yungel 



son et al. 1993). It is more likely that a less massive white 



WD undergoes collaps e due to electron captur e without 
exploding as an SN la (Nomoto & Kondo 1991). 

To examine the effec t of the stripping parameter, rjcff, 



we show, in Figure A13, the regions of SN la explosion for 
77eff = 0.3. Here, four cases of Mwd.o = 0.9, 1.0, 1.1, and 
I.2M0 are depicted. The region of Mwd,o = l.OAf0 for 
VcS = 0.3 is one third, in area, of the region for r]cS = 1 
see the region enclosed by the dash-dotted line in Fig. 
A12). In the limiting case of 77off = 0, we have again the 
constraint of g < 1.15 as in HKN96's modeling. 

4. DISCUSSION 

4.1. Recurrent novae as progenitors of SNe la 

First, we introduce a few binary systems which arc well 
understood in relation to our SN la progenitor model. T 
Coronae Borealis (T CrB) and RS Ophiuchi (RS Oph) are 
recurrent novae, which are binaries consisting of a very 
massive white dwarf a nd a lobe-filling re d-giant wit h or 



dwarf accretes hydrogen-rich matter from a red-giant com- 
panion and grows up to near the Chandrasekhar limit. If 
we include the mass-stripping effect by the strong WD 
wind, we easily reproduce the present states of T CrB and 
RS Oph systems. 

Specifying the initial parameters of AfwD,o — I.OM0, 
A^RG.o — I.3M0, and Pq — 135 d, we obtain the present 
state' of Mwd = 1.37M0, Mrg = O.7IM0 (with a he- 
hum core of Mhc = O.35Af0), P = 228 d, and the mass 
transfer rate of M2 '^ 1 x lO^^Af© yr~^. This set of the 



bital periods of 228 d (|Lines et al. 1988| ) and 460 d ( |Do 
brzycka fc Kenyon 1994| ), respectively. Two interpreta- 



initial parameters is shown in Figure All by a star mark 
(•) . It seems that T CrB is a critical (failing/succeeding) 
system for SN la explosion because this initial model 
corresponds to the lower boundary of the SN la region 
for A/wD.o = l.OA^0 (Fig. ^l|). It should be noted 
here that a recent analysis by Belczyhski & Mikolajew- 
ska (1998) shows, con trary to the previous results (e.g., 

the mass ratio of T CrB to be 
which implies a low-mass binary 



Webbink et al. 
A^rg/A/wd 



1987) 



tions have been proposed so far on the binary nature of T 
CrB: one is an episodic mass transfer event (model) onto a 



main-sequence s tar from a red-giant companion (e.g. Web- 



bink et al. 1987) and the other is a thermonuclear runaway 
event (model) on a mass-accreting wh ite dwarf as massiv e 
as the Chandrasekhar mass limit (e.g., ^elvelli et al. 199^ . 
Ultra-violet lines in its quiescent phase observed by lUE 
strongly indicate the existence of a mass- accreting whit e 
dwarf instead of a main-sequence star ( ^elvelli et al. 199^ ). 
Very rapid decline rates of their light curves indicate 
a very massive white dwarf close to the Chandrasekhar 
fimit , that is, A%d ^ 1.37 - 1.38Af0 for T CrB (|Katc 



q = AiRG/A^WD ^ 0.6 

system and is very consistent with the present numerical 

results. 

For RS Oph, if we start the calculation with A/vvd,o = 
l.OAf0, A^RG.o = 1.15A^0, and Pq = 240 d, we obtain the 
present state' of A^wd = 1.36A^0, ATrg = O.GOMq (with 
a helium core of Afno = O.39A^0), P = 460 d, and the 
mass transfer rate of M2 ^ 1 x lO~^Af0 yr~^. RS Oph 
also seems to be a critical system for SN la explosion. It 
should be noted that these sets of the initial parameters 
are not unique. 

U Scorpii (U Sco) is also one of the well-known recur- 
rent novae. Light curve fitting indicat es a very m assive 
white dwarf of AfwD = 1.37 - 1.38A'f0 (iKato 199C|. 1995 
1999). The orbital period is P = 1.23 d (pchaeferfc Ring- 



calculated the nova hght curves tor the white dwart masses 
of 1.2, 1.3, 1.35 and 1.377 Mq, where 1.377 Mq was cho- 
sen as a limiting mass being just below the mass at the 
SN la ex plosion in the W7 model (1.378Af0, [Nomoto et 
al. 1984|). Kato found that the light curve of the 1.377M0 



1995|,|1999). It should be noted that Kato (1995, 1999) wald 1995|), thus corresponding to the WD+MS system. 



Observations have suggested th at the companion of U Sco 



model is in better agreement with observational light curve 
of T CrB than the lower mass models. 

Very recently, other observational supports for a mas- 
sive white dwarf in T CrB have been reported: one is 



is extremely helium-rich (e.g., Williams et al. 1981 ), al- 
thou gh its compan i on is a slightly evolved ma in-sequence 
star ( ^chaefer 199C ; lohnston fc Kulkarni 1992 ). Its evolu- 
tionary path has long been regarde d as a puzzle in the the - 
ory of close binary evolution (e.g., Webbink et al. 1987 ). 
Very recently, Hachisu fc Kato (1999b) have elucidated the 
reason why its companion has a helium-rich envelope and. 
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at the same time, why the white dwarf is so massive as the 
Chandrasekhar mass hmit. 

In Hachisu & Kato's (1999b) U Sec scenario, they start 
from a progenitor binary system of ~ 7 — 8Mq and ~ 2Mq 
stars with the initial separation of ^ 500 R(^. The pri- 



mary component has first evolved to fill its Roche lobe 
when the helium core grows to ~ 1.4 — 1.6AIq. The bi- 
nary undergoes a common envelope evolution and shrinks 
to the separation of ~ lOi?© between the naked helium 
core and the '^ 2Mq main-sequence star. The helium core 
evolves to fill its Roche lobe and stably transfers almost 
pure helium onto the secondary because of the mass ra- 
tio M1/M2 ^ 0.79. As a result, the secondary becomes a 
helium-rich star. 

After the helium envelope of the primary is exhausted, 
the primary becomes a carbon-oxygen (C-l-0) white dwarf 
of 0.9 — I.OM0 and the secondary grows in mass to 
^ 2.5Mq. When the secondary slightly evolves to fill 
its Roche lobe, it transfers helium-rich matter onto the 
C+0 white dwarf on a thermal time scale. The white 
dwarf burns hydrogen atop the surface at a critical rate of 
Mcr ^ 2.0 X 1O-^(Mwd/M0 - 0.40) Mq yr^^ for helium- 
rich matter and blows excess matter in winds. The white 
dwarf now grows to near the Chandrasekhar mass limit 
and the mass transfer rate decreases to a few to several 
times 1O~^M0 yr~^. These pictures seems to be very con- 
sistent with the present distinct observational aspects of 
USco. 

To summarize, the recurrent novae seem to be a crit- 
ical system for SN la explosion. Recurrent novae are 
morphologically divided into three groups; dwarf compan- 
ion, slightly evolved main-sequence companion, and red- 
giant companion (e.g., Schaefer fc Ringwald 1995 ). T CrB 
(Porb == 228 d) and RS Oph (Porb = 460 d) belong to the 
last group of red-giant com panion. U Sco (Porb = 1-23 d; 
Schaefer fc Ringwald 1995|) and V394 CrA (Po^b = 0.758 d; 



Schaefer 199C ) belong to the middle group of the slightly 



evolved main-sequence companions. Two of the three sub- 
groups of the recurrent novae correspond to our progen- 
itors (WD+MS/WD+RG systems). This close relation 
between the recurrent novae and our progenitors strongly 
support our scenario of SN la progenitors. 

4.2. Yungelson & Livio's criticism 

Based on the population synthesis analysis, Yungelson fc 
Livio (1998) claimed that almost no realization frequency 
is derived for the original HKN96's WD-I-RG model. First, 
we briefly explain their analysis why the original model by 
HKN96 does not produce enough number of SNe la. Sec- 
ond, we point out that very wide binaries with the initial 
separation of a^ ^ 1500 Pq, which were not included in 
Yungelson fc Livio's (1998) analysis, are essentially impor- 
tant in our SN la modeling. 

A more massive component (mass of Mi,i) of a binary 
first evolves to a red-giant (AGB stage) and fills its in- 
ner critical Roche lobe. After a common envelope phase, 
the more massive component leaves a C-l-0 WD and the 
separation of the binary decreases by a factor of 



a. 



acE 



/ A%d \ 



M2 



V Ml,, J \Mi,,^MwD 



(36) 



evolutions, a/ (a;) the final (initial) separation, and M2 
the mass of the secondary. Adopting a standard value of 
"CE = 1, we obtain Uf/ai ~ 1/40—1/50 for Mwd ~ 1 Mq 
and M2 ~ 1 Mq, because a '^ 1 Mq WD descends from 



a main-sequence star of Mi,. 



mann 1986; Yungelson et al. 



7 

1995| ). 



8 Mq (e.g., Weide- 



Yungelson & Livio 

(1998) assumed that the separation of interacting binaries 
is smaller than Ui < 1500 Pq. Then, the most wide bi- 
naries has the separation of a/ ^ 30 — 40 Rq after the 
common envelope evolution. Its orbital period is Pq <j 20 
d for MwD,o ~ 1 Mq and Mrg.o ~ 1 Mq. There is no SN 
la region of the W D-I-RG systems for Pq ^ 20 d as seen 



from Figure All. Thus, they concluded that we cannot 
expect any SN la explosions from the right SN la region 
(WD+RG system) of Figure [All . 

If the WD+RG evolution starts only from an initial con- 
dition of Pq ^ 20 d and Mwd,o ~ IMq, however, the 
present states of T CrB or RS Oph cannot be reached be- 
cause of low mass transfer rates of |Mt| < 1 x 10~'^Mfn 



yr~^ (see also Figs. All and [A12| ). Thus, the existence of 
recurrent novae T CrB and RS Oph seems to be against 
the above Yungelson fc Livio's conjecture. The reason 
why T CrB or RS Oph are failed to be reproduced in Yun- 
gelson fc Livio's modeling is due to their assumption of 
ai <^ 1500 Rq. In what follows, we show that such wide 
WD+RG binaries as Pq '- 100 - 1000 d are born from 
initially very wide binaries with ai ~ 1500 — 40000 Rq. 

A star with the zero-age main-sequence mass of Mi^i <^ 
8 Mq ends up its life by ejecting its envelope in a wind of 
relatively slow velocities {v ^ 10 — 40 km s^"'^). These wind 
velocity is as low as the orbital velocity of binaries with 
the separation of Oi ~ 1500 - 40000 Rq for Mi,i ~ 7 Mq 
and M2^i ~ 1 M©. When the wind velocity is as low as the 
orbital velocity, the numerical factor i^ in equation (11) 
increases to 



1.7-0.55 



aflr. 



(37) 



mainly because outfiowing matter can get angular momen- 
tum from the binary motion by torque during its journey 
(see Appendix A). Here, v is the radial component of the 
wind velocity near the inner critical Roche lobe and the 
limiting case of ^w = 1-7 for w = was obtained by Nariai 
(1975) and Nariai fc Sugimoto (1976) for a test particle 
simulation ejecting from the outer Lagrangian points and 
Sawada, Hachisu, fc Matsuda (1984) for a 2-D (equato- 
rial plane) hydrodynamical simulation blowing a very slow 
wind from the primary surface which fills the inner critical 
Roche lobe. 

Combining the two expressions, we obtain 



1.7-0.55 



aQ 



orb 



(38) 



which is a good approximation to £^ in the region of the 
wind velocity from zero to infinity. Switching from equa- 
tion (n2) to equation ( |37| ) occurs at 1; ~ 1.5ariorb for q = 2. 
If we assume that slow winds blow from the primary and 
the mass of the secondary does not change (M2 = 0), the 
separation is calculated from 



where acE is the efficiency factor of common envelope 



a _ Ml + M2 
a ~ Ml + M2 



2^-2^ + 2:^ 
Ml M2 J 



Hachisu, Kato, and Nomoto 
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2L 



M1 + M2 



Ml 



M, 



Afi 



Ml 
M^ 



Ml 
'M2 



Ml 
Ml' 



M2 

(Ml > M2) 



(39) 



Therefore, the critical value of £^ is about 0.5q — 
0.5(M2/Mi) for shrinking/expanding of the separation. 
When V '^ afloih , the systemic loss of the angular momen- 
tum is estimated as i?w ~ 1-7 — 0.55{v/aQorh)^ '^1.0 from 
equation (||). We therefore have a/ a » 2M1/A/2 < 0. 

Once the binary system begins to shrink, its evolu- 
tion becomes similar to a common envelope evolution (see 
stages B-D in Fig. Al). As the separation shrinks, the 
orbital velocity of afJorb increases. If the wind velocity 
is almost constant, the ratio of ti/afJorb in equation (38) 
becomes smaller and smaller and the shrinking is accel- 
erated more and more. Thus, the separation is reduced 
by a factor of 1/40 - 1/50, i.e., a/ - 30 - 1000 Rq 
for Mi,i ~ 7 Mq and M2,i '- 1 Mq. The orbital pe- 
riod becomes Pq ~ 15 — 3000 d for AfwD,o '^ 1 Mq and 
M2 ~ 1 Mq {stage D in Fig. ^l]). These initial sets of the 
parameters are very consistent with the initial conditions 
of our WD-f RG progenitor systems. 

To summarize, we must include binaries with the sepa- 
ration of 



Qi < 5500R, 



■o 



1.7 



Ml, 



M2 



Mr 



10 km s-^ 



into the category of interacting binaries, where v is the 
velocity of slow wind (super wind) at the end of AGB evo- 
lution and ^^1.5 — 1.7 is a numerical factor defined by 



o q(l + 2q) 

^2 = 3.0-0.9-^ ^' 



(l + # 



(41) 



Here, we take the critical wind velocity of shrink- 
ing/expanding as 

Wcr = ^aflorh, (42) 

at the beginning of super wind phase. 

4.3. SN la frequency 

We estimate the SN la rate in our Galaxy coming from 
our WD-I-RG/WD-I-MS systems by using equation (1) of 
Iben & Tutukov (1984), i.e.. 



iy = 0.2- Aq- 



Mb 



Ma 



dM 



AlogA yr \ 



(43) 



where Aq, A log A, Ma, and Mb are the appropriate 
ranges of the mass ratio, of the initial separation, and the 
lower and the upper limits of the primary mass for SN la 
explosions in solar mass units, respectively. The estimated 
rate of the WD+RG/WD+MS systems is clos e to the ob- 
served rate in our Ga laxy, ly r^ 0.003 yr"-'^ (e.g. , Cappellarc 
et al. 1997; see also Yungelson fc Livio 1998), as will be 



shown below. 

4.4. WD+RG systems 

For the WD+RG progenitors, we assume that the ini- 
tial region of the separation includes a^ ~ 1500 — 40000 Rq 



as well as at < 1500 Rq (see discussions above in §4.2). 
Dividing the initial white dwarf mass of Mwd.o into four 
intervals, i.e., O.8-O.9M0, O.9-1.OM0, l.O-l.lM©, and 
1.1 — 1.2Mq, we estimate the realization frequencies. The 
mass range of AfvvD,o < O.7M0 is not included because 
no SN la explosions are expected for such low initial mass 
WDs. We omit the range of Mwd,o = O.7-O.8M0 because 
its realization frequency is too s mall to c ontribute to the 
SN la rate as seen in Figures Al2| and [Al3 . To estimate the 
initial separation(or orbital period), log A, and the initial 
lower/upper masses. Ma and Mb, of our WD-fRG sys- 
tems, we need to obtain the zero-age main sequence mass 
of the primary component [Mi^i) and the contraction fac- 
tor after the first common envelope phase. In single star 
evolutions, 0.7 — 1.2Af0 white dwarfs descend from stars 
with the zero-age main sequence mass of Mi '^ 3 — 8M0, 
i.e.. Ma ~ 3M0 and Mb ~ 8M0. More precisely, using 
Yungelson et al.'s (1995) equation (11) gives the final core 
mass (C-l-0 WD mass) vs. the zero-age main sequence 
mass relation, 



log- 



A/, 



c+o 



M 



-0.22 + 0.36 log 







Mo 

Mr, 



2.5 



(44) 



as numerically summarized in Table A4 , where Mwd,o = 
Afc-i-o ■ The initial separation should be larger than 

in order that the C+O core grows up to Afc+o — Mwd,o- 
Here f{q) « 0.5 for q = Mi^i/M2,i ~ 2 - 7 and the radius 
of stars at the AGB phase is given by 



^ii^^ ^ 1050 f^^^- 0.5 

Rq \ Mq 



0.68 



(46) 



(Iben & Tutukov 1984). For example, the initial separation 
should be larger than a^ ^ 1, 200 Rq for Mwd,o = 1-OMq 



as summarized in Table A4 

For the binary of AfwD,o = I.OM0 and M2,i = IMq, 
the contraction factor is estimated to be 1/37 by assum- 
ing the common envelope efficiency factor of acE ~ 1- 
The range of the separation after common envelope evo- 
lution becomes a/ ~ 32 — 1,120 Rq (corresponding to 
Po ~ 15 - 3, 070 d) because of a, - 1, 200 - 41, 700 Rq 
as summarized in Table K4. The orbital period of Pq ~ 
15 — 3, 07 d c overs the SN la region (WD+RG system) 
of Figure All . The bina ry p arameters for other Mwd,o 
are summarized in Table |A4| . The regions of the orbital 
period, logPo, also c overs me SN la region (WD+RG sys- 
tem) of Figure A12. Here we assume w = 10 km s^^ and 
V < VcT — C'^i^orb to meet the condition for the binary to 
contract (§4.2). Since the region of log A (logai) is shifted 
in parallel to the region of log af by the contraction fac- 
tor, the probability frequency for A log A is the same as 
for A log af . Then we approximately set as 



A log A « A log Pq 



3' 



(47) 



whe re A logPo is taken from the SN la region in Fig- 
ure A12 and the factor of 2/3 comes from the conver- 
sion between the period and the separation. Substitut- 
ing A log A = 0.6 • 2/3, Aq = 2.6/4.48 - 1.2/5.60 = 0.37, 
Ma = 4.48, Mb = 5.60, we obtain javd,o.8-o.9 = 0.0006. 
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The SNe la rates for other WD mass intervals are sum- 
marized in Table ^^. Then, the summation of SN la 
rates for three intervals (0.8 - 0.9Mq, 0.9 - I.OMq, and 
1.0 — I.IMq) gives z^RG — 0.0017 yr^^, which is large 
enough to explain the dominant part of the SN la rate 
in our Galaxy. If we further include the WD mass r ange 
of Mv^D.o = 1-1 ^ 1-2Mq, which is not shown in Table A5, 
the realization frequency increases to i/j^q = 0.0022 yr^^. 
Here, we have not shown the region for Mwd,o = 1.2Mq 



4) 



in Figure A12 because Webbink et al.'s (1983) empirical 
formula is valid for M2,o < 2.5 — S.OMq with a degenerate 
helium core. The range of Af-yvD.o — 1.2M0 exceeds this 
limit. 

To examine the effect of the stripping parameter, rj^^ii, 



After helium exhaustion, a carbon-oxygen core de- 
velops. When the core mass reaches 0.9 — I.OMq, 
the helium star evolves to a red-giant and fills again 
its inner critical Roche lobe. Almost pure helium is 
transferred to the secondary because the mass trans- 
fer is stable for the mass ratio q — M1/M2 < 0.79. 
The mass transfer rate is as large as ~ 1 x IO^^Mq 
yr"""^ but the mass-receiving main-sequence star (~ 
2 — 3M(T)) does not expand for such a lo w rate (e.g., 



Kippenhahn & Meyer-Hofmeister 1977) 



we est u2^^^^^u^^ji^^^^^^^^^^^^££^^^^^^ 



r*\r tm- 



ih(^ 



n "i 



oc gnm- 



5) The secondary has received 0.1 — OAMq (almost) 
pure helium and, as a result, it bec omes a helium- 
rich s tar as observed in U Sco (e.g., Williams et al, 
19811; iBarlow et al. 1981| ; |Hanes 1985^ [Sekiguchi et 



mar\7h.r\ in Table A£ The parameter region shrinks to one 



third in area compared with the case of rjcs = 1, so that 
the realization frequency is reduced to about one third of 
1 case, i.e., ly^G — 0.0008 yr^^. 



VeS 



4.5. WD+MS systems 



For the WD+MS progenitors, HKNU98 have found a 
new evolutionary path, which has not been taken into ac 



count in the previous works (e.g ., Rappaport et al. 1994 ; 



Di Stefano fc Rappaport 1994 ; Yungclson et al. 1996 



Yungelson & Livio 1998), and estimated the realization 
frequency to be as large as j^ms = 0.001 yr~^. We briefly 
follow their new evolutionary path in the following and dis- 
cuss the total rate of SN la explosions: If a ~ IMq C-I-0 
WD is descending from an AGB star, its zero- age main 
sequence mass is ^ 7Mq by equation ( [44[ ) and the binary 
separation is larger than at ~ 1350 Rq if the secondary 
mass is ~ 2Mq. Its separation shrinks to a/ ^ 70 Rq 
after the common envelope evolution with acE = 1 (the 
secondary mass of ~ 2Mq). Then the orbital period be- 
comes Pq ^ 4 d and too long to become an SN la as 
seen in Figure All. Therefore, the C-I-0 WD comes not 
from an AGB star having the radius of equation ( p^ ) but 
from a helium star whose hydrogen-rich envelope has been 
stripped away in the flrst common envelope evolution (at 
the red-giant phase with a helium core). Then, we follow 
the evolution of a binary consisting of a helium star and a 
main-sequence star. 

To estimate the decrease in the separation after the com- 
mon envelope phase, we use the radius to helium core mass 
(i?i — AfiHo) relation, which are taken from tables given 
by Bressan et al. (1993). As an example, let us con- 
sider a pair of 7Mq + 2.5Mq with the initial separation 
of Oi ~ 50 — 600i?Q . The binary evolves to SN la through 
the following stages: 

1) When the mass of the helium core grows to I.OMq < 
Ml. He < 1.4M0, the primary fills its Roche lobe and 
the binary undergoes a common envelope evolution. 

2) After the common envelope evolution, the system 
consists of a helium star and a main-sequence star 
with a relatively compact separation of a/ ~ 3 — 
40 Rq and Porb - 0.4 - 20 d. 

3) The helium star contracts and ignites central helium 
burning to become a helium main-sequence star. 
The primary stays at the helium m ain-sequence for 
~ 1 X 10'^ yr (e.g., [Paczynski 1971]). 



al. 1988|). The separation and thus the orbital period 



gradually increases during the mass transfer phase. 
The final orbital period becomes Porb ^ 0.5 — 40 d. 

6) An SN la explosion occurs when Pq = 0.4 — 5 d and 
M2 = Mms.o -- 2-3M0 in the logPo-Md,o plane for 
the syste m of Afi = Mwd,o ^ 0.9 — I.OM0 as seen 
in Figure A12. 



Therefore, the above pair of 7M0 -I- 2.5M0 can be a 
progenitor of SNe la if the initial separation is between 
50 — 15Oi?0, which initiates a common envelope evolution 
at the helium core mass of Mihc — 1-0 — I.2M0, corre- 
sponding to the initial orbital perio d of P prb.o = 0.5 — 5 d 
for the WD-I-MS systems in Figure All[ . In this case, we 
have A log A = log 150 - log 50 = 0.5 and q = 2.5/7 = 0.36 
in equation (p3|). 

Calculating 25 pairs of AIi,i and A/2,i, IIKNU99 have 
obtained the parameter region of SN la explosion as Aq = 
0.4, A\ogA = 0.5, Ma = 5.5^nd Mb = 8.5. Substituting 
these values into equation (43), we obtain the SN la rate 
of j/MS = 0.001 yr-i. Thus HKNU99 have shown that the 
frequency of the WD-I-MS systems is about one third of 
the inferred rate in our Galaxy, which is much larger than 
that of Yungelson & Livio's (1998) estimation. It should 
be noted that Yungelson & Livio (1998) have obtained the 
birth rate of 1 x 10~^ yr~^ for their models 15 and 16 
by relaxing all the constraints on the mass ratio of their 
binary models, although it is not a realistic case. 

The orbital velocity of the WD-I-MS systems is much 
faster than that for the WD-I-RG systems, i.e., ailorb ^ 
400 km s~i for Mwd = I.OM0 and Mms = 2.0Mq at the 
zero- age main sequence. Then, the switching from equa- 
tion (12) to equation ( |37| ) occurs at i; ~ 1.5ariorb ^ 600 
km s^ . This means that the wind velocity has to be faster 
than ~ 600 km s^^ in order to avoid the formation of a 
common envelope. Otherwise, winds carry large specific 
angular momentum and drastically shorten the separation 
to enhance the mass transfer and to eventually form a com- 
mon envelope. Our strong winds satisfy the condition of 
V ^ 600 km s~^ and this supports Li & van den Heuvel's 
evolutionary process and also HKNU99's modeling on the 
WD-I-MS systems. 

To summarize, the contribution of the WD+MS sys- 
tems to the SN la rate in our Galaxy is about one third 
of the inferred rate in our Galaxy. The total SN la rate 
of the WD+MS/WD+RG systems becomes t'Rc+MS = 
f^RG + f^MS = 0.003 yr^^, which is close enough to the 
inferred rate of our Galaxy. 
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4.6. Observational detection of hydrogen 

In our scenario, the WD winds form a circunistellar en- 
velope around the binary systems prior to the explosion, 
which may emit X-rays, radio, and Ha lines by shock heat- 
ing when the ejecta collide with the circumstellar envelope. 
In HKN96's model of binary evolution, the mass accretion 
rate decreases well below 1 x 10~^ Mq yr~^ as the 
white dwarf mass gets close to the Chandrasekhar limit. 
Thus the strong wind has ceased when the white dwarf ex- 
plodes (HKN96). In contrast, the mass accretion rate in 
the present models is still as high as 1 x 10~^ Mq yr"-'^ for 
some of the white dwarfs near the Chandrasekhar limit, so 
that such a white dwarf explode in the strong wind phase. 

Our strong wind model of case PI predicts the presence 
of circumstellar matter around the exploding white dwarf. 
Whether such a circumstellar matter is observable depends 
on its density. The wind mass loss rate from the white 
dwarf near the Chandrasekhar limit is as high as Af ~ 1 
X 10~^ —1 X 10~^ MfTi yr~^ and the wind velocity is u = 
1000 km s^^ (Fig- |A5| ). Despite the relatively high mass 
loss rate, the circumstellar density is not so high because 
of the high wind velocity. For steady wind, the density is 
expressed by M/v (= iwr'^p). Normalized by the typical 
red-giant wind velocity of 10 km s~^, the density measure 
of our white dwarf wind is given as M/wio '^ 1 x 10~^° 
— Ix 10^^ Mq yr~^, where wio = u/10 km s^^. 

Behind the red-giant, matter stripped from the red-giant 
component forms a much dense circumstellar tail. Its rate 
is as large as ^ 1 x lO^^Af0 yr^^ with the velocity of ~ 100 
km s""'^. The density measure of the dense red-giant wind 
thus formed is M/vio ^ I x 10"^ Mq yr'^ 

Further out, circumstellar matter is produced from the 
wind from the red-giant companion, which is too far away 
to cause significant circumstellar interaction. 

For cases P2 and P3, winds from the WD have stopped 
before the explosion. Therefore, circumstellar matter is 
dominated by the wind from the red-giant companion 
whose velocity is as low as ^ 10 km s~^. 

At SN la explosion, ejecta would collide with the cir- 
cumstellar matter, which produces shock waves propagat- 
ing both outward and inward. At the shock front, particle 
accelerations take place to cause radio emissions. Hot plas- 
mas in the shocked materials emit thermal X-rays. The 
circumstellar matter ahead of the shock is ioni zed by X- 
rays and pro duce recombination Ha emissions ( [Cumniing 
et al. 1996 ). Such an interaction has been observed in 



Type lb, Ic, and II supernovae , most typically in SN 1993 J 



e.g., [Suzuki fc Nomoto 1995| and references therein) 



For SNe la, several attempts have been made to detect 
the above signature of circumstellar matter. There has 
been no radio and X-ray detections so far. The upper limit 
set by X-ray observations of SN 1992A is M/viq = (2 — 3) 
X 10-^5 Mq yr-i ( |Schlegel fc Petre 1993|) . Radio obser- 
vations of SN 1986G have provided the most strict upper 
limit to th e circumstellar d ensity as M/vio = 1 x 10^^ 
Mq yr-1 ( |Eck et al. 19951 ). This is still 10 - 100 times 
higher than the density predicted for the white dwarf wind 
for case PI. For cases P2 and P3, if the wind mass loss 
rate from the red-giant is significantly higher than 1 x 
10^^ Mq yr~^, radio detection could be possible for very 
nearby SNe la as close as SN 1986G. (Note also that SN 
1986G is not a typical SN la but a subluminous SN la.) 



For Ha emissions, Branch et al. (1983) noted a small, 
narrow emission feature at the rest wavelength of Ha, 
which is blueshifted by 1800 km s~^ from the local inter- 
stellar Ca II absorption. Though this feature was not ob- 
served 5 days later, such high velocity hydrogen is expected 
from the white dwarf wind model. For SN 1994D, Cum- 
ming et al. (1996) obtained the upper hmit of M/vio — 6 
X 10~^ Mq yr~^. Further attempts to detect Ha emis - 
sions are highly encouraged (Lundqvist & Cumming 1997). 



5. CONCLUSIONS 

Progenitors of SNe la have not been identified yet either 
theoretically or observationally. In the present paper, we 
propose a new evolutionary process which produces a wide 
enough symbiotic channel to SNe la. In this channel, the 
white dwarf accreting mass from the lobe-filling evolved 
companion (red-giant) grows to the Chandrasekhar limit 
and explodes as an SN la. In what follows we summarize 
our main results: 

1. In HKN96, we have found the most important mech- 
anism in our symbiotic channel to SNe la: when the 
mass accretion rate onto the white dwarf exceeds 
a critical rate of Mcr « 0.75 x 1O"^(Mwd/M0 - 
0.40)Mo yr~^, a strong optically thick wind blows 
from the white dwarf. It stabilizes the mass transfer 
even if the red-giant has a deep convective envelope. 
HKN96 have shown that the mass transfer becomes 
stable only when the mass ratio of RG/WD is smaller 
than 1.15. In the present study, we have found that 
this constraint is removed when the effect of mass- 
stripping from the red-giant envelope by the strong 
wind is large enough. Thus the symbiotic channel 
produces SNe la for a much (more than ten times) 
wider range of the binary parameters than that of 
HKN96's estimation, thus being able to account for 
the dominant part of the inferred rate of SNe la in 
our Galaxy. It should be noticed, however, that the 
realization frequency depends on the efficiency factor 
of mass-stripping T^off and the realization frequency 
is consistent with the inferred rate in our Galaxy for 
77eff ~ 0.3 - 1.0. 

2. Yungelson fc Livio (1998) estimated the realization 
frequency of HKN96's original model and concluded 
that almost no realization frequency is derived be- 
cause the binary shrinks too much (Porb < 20 d) to 
produce an SN la after the first common envelope 
evolution. In their population synthesis code, how- 
ever, Yungelson & Livio (1998) assumed that the 
initial separation is smaller than ~ 1500 Rq, i.e., 
ai ^ 1500 Rq. This assumption neglects the effect 
of slow winds at the end of stellar evolution. When 
the velocity of the slow wind is as small as the or- 
bital velocity, the slow wind gets angular momentum 
through torque by the binary motion and, as a result, 
the binary shrinks so much. When the wind velocity 
is as slow as '^ 10 km s~^, we must include the bina- 
ries with the separation of a^ ~ 1, 500 — 40, 000 Rq. 
Such an extremely wide binary will shrink to an ap- 
propriate range of the separation (20 d < Pq < 3000 
d for a pair of Mwd.o ~ 1-^^© and M2,i ~ l-A^o) 
after the slow wind (super wind) and the common 
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envelope evolution. If this effect is included in the 
calculation of SN la rate, we have a reasonable value 
of the realization frequency, i^rg = 0.002 yr""'^, for 
our white dwarf plus red-giant (WD+RG) systems. 

3. Yungelson & Livio (1998) estimated the realization 
frequency of Li & van den Heuvel's (1997) white 
dwarf plus main sequence star (WD+MS) model and 
concluded that the total frequency of their mod- 
ified HKN96's model and Li & van den Heuvel's 
model does not exceed 0.0002 yr~^, a tenth of the 
inferred rate. However, we believe that an impor- 
tant evolutionary path was not included in Yungel- 
son & Livio's (1998) analysis, that is, the primary's 
helium star phase: the primary becomes a naked 
helium star after the first common envelope phase 
if the mass transfer begins at the phase of a red- 
giant with a helium core. This helium star eventu- 
ally leaves a C-FO WD by transferring the helium 
envelope to the secondary. As a result, the sec- 
ondary may have an ext remely helium-rich envelope 



such as in U Sco (e.g., Williams ct al. 1981 



low ct al. 1981 



Hanes 1985 



Sckiguchi ct al. 



Bar 

1988^ 



'I'liis evolutionary path indicates that the secondary 
of the WD-I-MS systems has a helium-enriched en- 
velope thus forming an accretion disk which shows 
strong helium lines as s een in the luminous super- 
soft X -ray sources (e.g., Kahabka fc van den Hcuvej 



1997 for a recent review). Including this evolutionary 
path in the estimation of SN la rate, HKNU99 have 
obtained a much larger frequency of z/ms = 0.001 
yr~^ for the WD-I-MS systems than that estimated 
by Yungelson & Livio (1998). Thus the total fre- 
quency of our WD+RG/WD4-MS systems is as large 
as J^ = t^RG + t^MS = 0.003 yr^^, which is consistent 
with the inferred rate in our Galaxy. 

4. Recurrent novae are morphologically divided into 
three groups; dwarf companion, slightly evolved 
main- sequence companion, and re d-giant companion 
(e.g., ^chaefer fc Ringwald 1995| ). T CrB (Porb = 
228 d) and RS Oph (Porb = 460 d) belong to the 
last group of red-giant companion. It has been ar- 
gued that their white dwarf masses are very close to 
the Chandrasekhar mass limit. These systems corre- 
spond to near the lower board er of our SN la region 
(WD+RG) as shown in Figure All and its evolution- 
ary path is reasonably understood by our WD+RG 



systems. On the other hand, U Sco (Porb — 1-23 d) 
and V394 CrA (Po^b = 0.758 d) belong to the mid- 
dle group of slightly evolved main-sequence compan- 
ion. Our WD+MS model naturally yields a helium- 
enriched envelope of the secondary star as well as a 
near Chandrasekhar mass limit white dwarf as has 
been observationally suggested. The evolutionary 
path of the WD+MS systems is also very consistent 
with the middle group of recurrent novae. 

The immediate progenitors in our symbiotic chan- 
nel to SNe la may be observed as a symbiotic star. 
The photospheric temperature of the mass-accreting 
white dwarf is kept around Tph ^ 1 x 10^ — 2 x 10^ 
K during the wind phase. The hot star may not be 
observed in X-rays during the strong wind phase due 
to self-absorption by the wind itself. Some progen- 
itors stop blowing the wind before the SN la explo- 
sion, thus being observed as a luminous supersoft 
X-ray source. In some progenitors, very weak hy- 
drogen shell flashes are triggered before the SN la 
explosion; such a progenitor may be observed as a 
recurrent nova like T CrB or RS Oph. 

Radio emission from the circumstellar gas is pre- 
dicted. If the white dwarfs explode in the strong 
wind phase, however, radio emission is lower than 
the current observational limit because the wind 
velocity is as fast as 1000 km s~^ and the den- 
sity of the circumstellar medium is too tenuous to 
be observed even if the wind mass loss rate is as 
large as ~ 1 x lO^^Af0 yr^^ or more. If the wind 
has stopped before the explosion, the circumstellar 
matter is dominated by the low velocity wind from 
the red-giant companion; then observations of radio 
emission would be easier. Detection of high velocity 
hydrogen feature from the strong wind is also pre- 
dicted. 
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APPENDIX 



To obtain the relation between u/afiorb and ^^ in equations (|1 
moving under the Roche potential and Coriolis force, i.e.. 



SPECIFIC ANGULAR MOMENTUM OF WINDS 

and (B9|) , we calculate many orbits of test particles 



(Px 


-1— 


1 X — Xi 


q x~X2 


df^ 


l + q rl 


l + q rl 


d'y 
di2 


^dx 


1 y 

l + qrl 


q y 

l + qrl' 


d'z 


1 z 


q z 
l + qrl' 




dt^ 


l + qrl 





(Al) 



Hachisu, Kato, and Nomoto 



13 



where {x, y, z) is the position of the particle, q — AI2/M1 the mass ratio, ri and r2 are the distances from the primary 
and from the secondary to the test particle, respectively, and calculated from 



2 / \2 

2 / \2 

r2 = (x- X2) 



(A2) 



Setting the center of gravity of the binary at the origin of the coordinates, we assume that the primary is located on the 
negative side of the x-axis {xi, 0, 0) and the secondary is located on the positive side of the x-axis {x2, 0, 0), i.e.. 



Xi 
X2 



q 



1 



1 



(A3) 



We also assume a = 1, G = 1, and r^orb = 1 (or Mi + M2 = 1) in our dimensionless form. 

The primary (AGB star) blows spherically symmetric winds with the initial velocity oivo- This process can be simulated 
by ejecting test particles from the primary surface at the radius which is much smaller than the inner critical Roche lobe, 
e.g., one tenth of the inner critical Roche lobe of the primary; i.e., we set 



vq, 



at n =0.1i?i. 



(A4) 



The trajectory of the wind may be approximated by a trajectory of the particle with the same initial velocity and position 
when the wind is supersonic and does not form a shock. Here, we assume the equatorial symmetry of the wind. Dividing 
the primary surface into 64 x 256 parts, i.e., the azimuthal angle (from (p — to (f> = 27r) into 256, {A(f> — 27r/256), and the 
inclination angle (from 6 — to 9 = 7r/2) into 64 {A0 — 7r/2/64), we eject test particles from each center of the surface 
elements with the initial radial velocity of vq. We attach the mass loss rate of vq sin 9 iA9 Ac/) /Air to each particle. 
The radial component of the wind velocity Vr near the inner critical Roche lobe is calculated from 



2E 



1 / dxi 



n 



dt 



{Xi 



dyi 
dt 



dzi \ sm9iA0A(j) 



dt 



in 



at ri^i = Rl, 



(A5) 



where the position of each test particle is denoted by {xi, Hi, Zi) and the radial means the direction from the center of the 
primary to the test particle. We estimate the average specific angular momentum of the test particles by 



e. 






dyi 
' dt 



dxi 



s\n9iA9A4) ^^-^ s\x\9iA9A(t) 



in 



-/2 






in 



at r, 



10. 



(A6) 



It should be noted that some test particles are trapped to the secondary and never reach the radius of r = 10 when the 
radial velocity is smaller than the orbital velocity, i.e., v^ < aSlqrh — 1. We do not include these particles in the calculation 
of the specific angular momentu m o f the wind in equation (A6). 

The integration of equation (Al) is based on the second order leap-frog method. Five cases of the mass ratio, q — 
M2/M1 = 3, 2, 1, 1/2, and 1/3, are calculated for various initial velocity vq. The relation between the specific angular 
momentum £w and the radial velocity near the inner critical Roche lobe of the primary v{= Vr) is plotted in Figure 



A14. When the radial velocity of the wind is faster than twice the orbital velocity, i.e., v > 2 (or v ^ 2af7oib), the 



specific angular momentum is approximated by the limiting value of equation (12) for very fast winds. For wind velocities 
lower than 2, the specific angular momentum rapidly increases. We find that the values are located approximately on a 
quadratic line of 

£w« l■7-0.55i;^ (A7) 

where the limiting case of £w — 1-7 for 1; = was obtained by Nariai (1975) and Nariai & Sugimoto (1976) for a test 
particle simulation and Sawada et al. (1984) for a two-dimensional (equatorial plane) hydrodynamic simulation. Thus, 
the specific angular momentum is approximated by 



1.7-0.55 



aQ 



orb 



1 



q 



(A8) 



at least for these five different mass ratios, although three-dimensional hydrodynamic calculations should be done in order 
to obtain a definite conclusion of the specific angular momentum of the winds. 
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Table Al 
Numerical factors of mass-stripping effect 



q 


0.5 


1 


2 


3 


5 


9{q) 


0.013 
0.144 


0.019 
0.025 


0.025 
0.006 


0.030 
0.044 


0.036 
0.105 
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Table A2 
Three typical cases of SN Ia progenitor evolution 



AfwD,o 

{Mq) 



Mrg.O 
{Mq) 



Po 
(day) 



Mt 
(Mo yr-i) 



(Mo yr-i) 



(Mq yr-l) 



PI 
P2 
P3 



1.0 
1.0 

1.0 



2.0 
1.6 
1.3 



300 
300 
300 



-8.7 X lO-'^ 
-6.2 X lO"'^ 
-5.2 X 10 



7 



-4.8 X lO-'^ 
-2.3 X lO"'^ 
-1.3 X 10 



7 



-2.0 X 10-6 
-1.0 X 10"" 
-5.4 X lO"'^ 



Table A3 
Three states of immediate SN Ia progenitors 



history 



SN la explosion 



PI WIND 

P2 WIND 
P3 WIND 



WIND 
WIND 

sss 



WIND 
SSS 
RN 



WIND 

SSS 
RN 



Table A4 
Contraction factor (acE = 1) 



-'W'wD.O 


Mi.o 


ai 


at 


Po 


M2 


contraction 


(Mq) 


(Mq) 


(Rq) 


(Re) 


(day) 


(Mq) 


factor 


0.7 


3.20 


730-21,900 


64-1,920 


46 - 7, 480 


1.0 


1/11 


0.8 


4.48 


910-29,300 


44-1,420 


25 - 4, 640 


1.0 


1/21 


0.9 


5.60 


1, 060 - 35, 800 


36-1,220 


18 - 3, 600 


1.0 


1/29 


1.0 


6.63 


1, 200 - 41, 700 


32-1,120 


15 - 3, 070 


1.0 


1/37 


1.1 


7.58 


1, 340 - 47, 200 


30-1,060 


13 - 2.740 


1.0 


1/45 


1.2 


8.48 


1, 460 - 52, 400 


28-1,020 


12 - 2.530 


1.0 


1/52 



Table A5 
Realization frequency of SNe Ia (?7off = 1) 

MwD A log A Ma Mb A<j i^wd 

(Mq) (Mq) (Mq) (yr-l) 



0.8-0.9 
0.9- 1.0 
1.0- 1.1 



0.6 
1.0 
1.5 



2/3 
2/3 
2/3 



4.48 
5.60 
6.63 



5.60 
6.63 

7.58 



0.37 
0.36 
0.36 



0.0006 
0.0006 
0.0005 



Table A6 
Realization frequency of SNe Ia (?7off' = 0.3) 



A^WD.O 

(Me) 


A log A 


Ma 
(Mq) 


Mb 

(Mq) 


M 


(y-1) 


0.9-1.0 

1.0- 1.1 

1.1- 1.2 


0.7-2/3 
1.1-2/3 
1.8-2/3 


5.60 
6.63 

7.58 


6.63 

7.58 
8.48 


0.21 
0.23 
0.23 


0.0002 
0.0003 
0.0003 
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Progenitor of Type la Supernovae 
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ro 



total (wind + nuclear) 



wind mass loss rate 



X=0.7 Z=0.02 
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-4 
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Fig. A2. — Wind mass loss rate Mwind (dashed line) and the total mass decreasing rate of the hydrogen-rich envelope on the white dwarf 
Mnuc + -'\^wind (solid line), i.e., nuclear burning rate Mnuc plus wind mass loss rate M^indi S'l'e plotted against the envelope mass AM for 
WDs with masses of O.GMq, O.SMq, I.OMq, 1.2Mq, 1.3Mq, and 1.377M0. White dwarf mass is attached to each hne. The metallicity is the 
solar value oi Z = Zq = 0.02. There exists only a static (no wind) solution below the break of each solid line while there exists only a wind 
solution above the break for a given envelope mass of AM. The optically thick winds blow when the mass transfer rate from the companion 
star IM2I = |Mnuc + M^indl exceeds Mcv = 0.75 X 10-^{Mwd/Mq - 0.4)Mq yr-l. 
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Progenitor of Type la Supernovae 
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Fig. A3. — Photospheric temperature Tpj, is plotted against the decreasing rate of the envelope mass Mnuc + M„i^^ for WDs with 
masses of O.GMq, O.SMq, I.OMq, 1.2Mq, I.SMq, and 1.377Mq. Same envelope models as in Fig. 2. Optically thick winds blow when 
logTph (K)<5.3. 
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'og Rph (R©) 



Fig. A4. — Photosphcric radius i?pij is plotted against the decreasing rate of the envelope mass Mnuc + Mwind for WDs with masses of 
O.6M0, O.8M0, I.OMq, 1.2Mq, 1.3Mq, and 1.377M0. Same envelope models as in Fig. 2. Optically thick winds blow when _Rph > O.IRq. 
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Fig. A5. — Photosphcric velocity Vp\i is plotted against the decreasing rate of the envelope mass Mnuc + Mwind for WDs with masses of 
O.6M0, O.8M0, 1.OM0, 1.2M0, l.SM©, and 1.377M0. Same models as in Fig. 2. 
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Fig. A6. — Ratio of the photospheric velocity to the escape velocity there iipij/vesc is plotted against the decreasing rate of the envelope 
mass Mnuc + M^ind for WDs with masses of O.GMq, 0.7Mq, O.SMq, 0.9Mq, I.OMq, 1.2Mq, 1.3Mq, and 1.377M0. Same models as in Fig. 
2, but 0.7Mq and 0.9Mq are added. We regard the wind as "strong" when the photospheric velocity exceeds the escape velocity there, i.e., 
t^ph > fesc- If not, it is regarded as 



We regard the wind as "strong 
"weak." 
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Progenitor of Type la Supernovae 




Fig. A7. — Fast winds from the white dwarf (WD) collide with the surface of the red-giant and strip mass from the red-giant. The normal 
component of the wind velocity to the red-giant surface is dissipated by forming a shock and it heats up the surface. Then, a part of the 
surface mass is ablated and blown in the wind. 
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Fig. A8. — Time evolution of an SN la progenitor system for case PI (explosion during the wind phase). The initial parameters arc shown 
at the top of the figure. The white dwarf mass increases to 1.38Mq and explodes as an SN la at i = 7.2 X 10^ yr. The solid lines show the 
masses of the white dwarf (Mwd) and the red-giant companion (Mrq)- The dashed lines show, from top to bottom, the net mass accretion 
rate onto the white dwarf, the wind mass loss rate, and the mass decreasing rate of the red-giant companion, respectively. 
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Fig. A9. — Same as Fig. 8 but for case P2 (explosion during the steady hydrogen shell burning phase). The strong wind stops at the time 
indicated by an arrow. The white dwarf mass increases to 1.38Mq and explodes as an SN la at t = 7.4 X 10^ yr. 
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Fig. AlO. — Same as Fig. 8 but for case P3 (explosion during the very weak shell flash phase). The strong wind stops at the time indicated 
by an arrow. The hydrogen shell burning becomes unstable to trigger very weak shell flashes at t = 6.5 X 10^ yr. The white dwarf mass 
increases to 1.38M0 and explodes as an SN la at t = 1.04 X 10^ yr. 
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Fig. All. — Final outcome of close binary evolution in the logPo — Mj q plane. Here, Mj q is the donor mass, Mms.o or Mrg.O- Final 
outcome is either an unstable mass transfer (H4{q) > 0) at the beginning (forming a common envelope; denoted by x), or an SN la explosion 
(denoted by ®, Qi or 0) or a nova (denoted by a open triangle), or a central helium flash (denoted by a filled triangle). ®: wind phase at 
SN la explosion (PI). O- wind stops before SN la explosion but the mass transfer rate is still high enough to keep steady hydrogen shell 
burning, i.e., |Mtj > Mat (P2). 0: wind stops before SN la explosion and the mass transfer rate is decreasing between Mlo^^, < jMt| < Met 
at SN la explosion (P3). The region producing an SN la is bounded by a solid line. The left/right region corresponds to the WD+MS 
(compact)/WD+RG (wide) system, respectively. A star mark (*) denotes an initial position of T CrB in the WD+RG systems. 
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Fig. A12. — The region to produce SNe la in the logPo - Md,o plane for five initial white dwarf masses of 0.75Mq, O.SMq, 0.9Mq, I.OMq 
(heavy solid line), and I.IMq. The region of MwD.o = 0.7Mq almost vanishes for both the WD+MS and WD+RG systems, and the region 
of MwD,o = O.75M0 vanishes for the WD+RG system. Here, we assume the stripping efficiency of rjcS = 1- For comparison, we show only 
the region of M^d o = I.OMq for a much lower efficiency of rj^g = 0.3 by a dash-dotted line. 
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Fig. A13. — Same as Fig. 12 but for the much lower mass-stripping efficiency of rjcS = 0.3. We add the region of MwD.o = 1.2Mq both 
for the WD+MS/WD+RG systems (dotted hnes). The region of Mwd.o = 0-SMq vanishes for the WD+RG system. 
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q=3.0 




v/aQ 



Fig. A14. — Specific angular irLomentum of wind is plotted against the outflowing velocity near the inner critical Roche lobe. The specific 
angular momentum and the outflowing velocity are measured in units of a^fl and af2, respectively, where a is the separation and £7 is the 
orbital angular velocity. Five cases of the mass ratio are examined, i.e., q = M2/M1 = 3, 2, 1, 1/2, and 1/3. It is assumed that the wind 
blows from the primary. The limiting case of £w = 1.7 for d = is taken from Nariai (1975), Nariai & Sugimoto (1976), and Sawada et al. 
(1984). 



